function [DWL_pertax,DWL_perinc]=get_DWL(estimates)

mu1=0;
mu2=NaN;
mu3=NaN;
sigma2=NaN;
sigma3=NaN;
fdist1='Normal';
fdist2='';
fdist3='';
alpha1=NaN;
alpha2=NaN;
alpha3=NaN;
K1=estimates.config.K1;
K2=estimates.config.K2;
K3=estimates.config.K3;

config=estimates.config;
p_lb=config.p_lb;
p_ub=config.p_ub;

y_lb=config.y_lb;
y_ub=config.y_ub;

q_lb=config.q_lb;
q_ub=config.q_ub;

y_interest_univ_DWL=config.y_interest_univ_DWL; 
p_interest_univ_DWL=config.p_interest_univ_DWL; 

q0=config.q0;

p_interest_DWL=kron(p_interest_univ_DWL,ones(size(y_interest_univ_DWL,1),1));
y_interest_DWL=kron(ones(size(p_interest_univ_DWL,1),1),y_interest_univ_DWL);
tau_vec=[0.25; 0.50; 0.75]; 

q_hat_est=berk_ml.invert_Ginv(estimates,q0,tau_vec,p_interest_DWL,y_interest_DWL);

show_txt=0;
flag_median=0;

[DWL_pertax,DWL_perinc]=show_DWL(p_interest_univ_DWL,y_interest_univ_DWL,p_interest_DWL,y_interest_DWL,tau_vec,K1,K2,K3,p_lb,p_ub,q_lb,q_ub,y_lb,y_ub,mu1,mu2,mu3,sigma2,sigma3,fdist1,fdist2,fdist3,alpha1,alpha2,alpha3,estimates,q_hat_est,show_txt,'','',flag_median);

%

